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STABLE MULTIVARIATE M^-EULERIAN POLYNOMIALS 

MIRKO VISONTAI AND NATHAN WILLIAMS 

Abstract. We prove a multivariate strengthening of Brenti's result that every root of the 
Eulerian polynomial of type B is real. Our proof combines a refinement of the descent 
statistic for signed permutations with the notion of real stability — a generalization of real- 
^— s , rootedness to polynomials in multiple variables. The key is that our refined multivariate 

^\i ' Eulerian polynomials satisfy a recurrence given by a stability-preserving linear operator. 

Our results extend naturally to colored permutations, and we also give stable general- 
C^ ' izations of recent real-rootedness results due to Dilks, Petersen, and Stembridge on affine 

Eulerian polynomials of types A and C. Finally, although we are not able to settle Brenti's 
real-rootedness conjecture for Eulerian polynomials of type I?, nor prove a companion con- 
jecture of Dilks, Petersen, and Stembridge for affine Eulerian polynomials of types B and 
D, we indicate some methods of attack and pose some related open problems. 

o 
u 

c^ ! 1. Introduction 

In this paper, we study the real-rootedness property of Eulerian polynomials for Coxeter 
groups from a combinatorial perspective. There is a well-known combinatorial interpretation 
of the Eulerian polynomial An{x) as the descent generating polynomial for permutations in 
the Coxeter group An, the group Sym(n-l-l) of all permutations on n + 1 letters (see p^ITT] ). 
fd ■ The notion of a descent can be extended to elements of all finite Coxeter groups as follows: 

O ■ for an element a of the Coxeter group W the descents are exactly those generators s oi W 

cn ! whose action on a reduces its length. In [7] , Brenti used this interpretation of a descent to 

define the VF-Eulerian polynomial, denoted W{x), for any finite Coxeter group W. 

Brenti showed that many classical results about An{x) hold for the other Eulerian poly- 
nomials as well. In this paper we will investigate the remarkable property that An{x) has 
only real roots, a result due to Frobenius p3]. Brenti proved the analogous result for type 
B, and checked by computer that it also held for the exceptional cases, but left type D — the 
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. ^.' only remaining case — as a conjecture. 

Conjecture 1.1 (Conjecture 5.2 in [7]). For every finite Coxeter group W, the descent 
generating polynomial W{x) has only real roots. 

Dilks, Petersen, and Stembridge later extended the definition of Eulerian polynomials to 
include affine descents, and proposed the following companion to Brenti's conjecture. 

Conjecture 1.2 (Conjecture 4.1 in [9]). For every finite Weyl group W , the affine descent 
generating polynomial W{x) has only real roots. 

Again, this was not completely proved — A„(x) and C„(x) were shown to have only real 
roots, the exceptional cases were verified, but the real-rootedness of the affine Eulerian 
polynomials of types B and D remains an open problem. 



In this paper, we build on the idea of real stability — a generalization of the notion of real- 
rootedness to multivariate polynomials. We combine this with simple recurrences for multi- 
variate refinements of certain Eulerian polynomials to provide simple proofs of multivariate 
generalizations of known real-rootedness results. Specifically, we give a general framework to 
show that the recurrence relations satisfied by multivariate W- and W^-Eulerian polynomials 
(for certain finite Coxeter groups W) are stability-preserving. We then use properties of 
stability to show that this implies that the univariate counterparts of these polynomials are 
also stable, which is equivalent the statement that they have only real roots. 

The remainder of this paper is structured as follows. In Section [21 we introduce notation, 
define the VF-Eulerian and the affine VF-Eulerian polynomials for finite Coxeter groups and 
finite Weyl groups, respectively. We also review the required definitions and results related 
to real stability. For clarity and completeness, we begin in Section 13.11 with a proof due to 
Branden of the stability of the multivariate Eulerian polynomial of type A. In Sections 13.21 
and l3.3l we generalize this idea in several directions simultaneously, to type B (signed permu- 
tations) and the generalized symmetric group (colored permutations), and also to multiple 
q variables. Sections 14.11 and 14.21 then address the affine Eulerian polynomials for types A 
and C. The unresolved cases of Conjecture 11.11 (type D) and Conjecture 11.21 (types B and 
D) are examined within our multivariate framework in Section [51 We conclude with a dis- 
cussion about the connection between our statistics with Catalan numbers, Motzkin paths 
and Laguerre polynomials. 



2. Preliminaries 

We begin by introducing some notation. For a positive integer n, let [n] be the set 
{1, . . . , n} and let x be the n-tuple Xi, . . . , x„; for example, x + y = xi + ?/i, . . . , x„ + ?/„. For 
7 a set (or multiset) with entries from [n], we let x^ = nieT-^JJ ^'^^ example, (x + y)'"' = 
nr=i('^« + yi). The cardinality of T is written |T|. We denote the concatenation of x and y 
by X, y or x; y. 

We apply a function f oi n variables to an n-tuple x by writing /(x) = f{xi, . . . , x„). We 
will often deal with functions that have x and y as variables, and so we define the special 
symbol d = Yl^=ii'£r ~^ W~^ ^^ ^ shorthand for the sum of partial derivatives with respect 
to all X and y variables. 

Finally, the theorems and propositions that are taken from previous works are clearly 
marked by a reference (indicating the source); as far as we know, all other results are new. 

2.1. ly-Eulerian Polynomials. Let S* be a set of Coxeter generators, m be a Coxeter 
matrix, and 

W={S: {ss'r^''''^ = e, for s, s' G S, m{s, s') < oo) 

be the corresponding Coxeter group (see [2]). Given such a Coxeter system {W,S) and 
a G W, we denote by iw{<^) the length of cr in VF with respect to 5*. 

Definition 2.1. For W a finite Coxeter group, the descent set oi a ^ W is 

'Dwia) = {seS: iw{(Ts) < ew{(^)}- 
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Definition 2.2. For W a finite Coxeter group, the W-Eulerian polynomial is the descent 
generating polynomial 

The above definitions were extended in [S] to include afiine descents. 
Definition 2.3. For W a finite Weyl group, the afftne descent set of cr G VF is 

i)wi(T) = DH^(a) U {so : iw{(^So) > iw{<7)}, 

where Sq is the refiection corresponding to the lowest root in the underlying crystallographic 
root system. See [9] for further details and the motivation behind this definition. 

Definition 2.4. For W a finite Weyl group, the W-Eulerian polynomial is the affine descent 
generating polynomial (over the corresponding finite Weyl group W) 

o-Giy 

2.2. Real Stable Polynomials. We define real stability, which generalizes the notion of 
real-rootedness from real univariate polynomials to real multivariate polynomials. 

Let IK+ = {z G C : Im(2;) > 0} denote the open upper complex half-plane and similarly 
let 'K^ = {zeC: Im(^) < 0}. 

Definition 2.5. A polynomial / G M[x] is (real) stable if / = or for any z G !K" , /(z) 7^ 0. 

Note that a univariate polynomial f{x) G M\x\ has real roots if and only if it is stable. 
Following [21], we let ©r[x] denote the set of stable polynomials in M[x]. 

In this paper, we have a fixed template for our proofs. We argue by induction, first 
checking stability (by hand) for the base case. Next, we establish recursive formulas of the 
following form 

where Wn is the multivariate VF-Eulerian polynomial of a group W of rank n, and T is some 
linear operator. Finally, we show that the linear operator T is stability-preserving, using the 
following theorem. 

Recall that a polynomial /(x) is multiaffine if the power of each indeterminate Xi is at 
most one. For a set 7 of polynomials, let [P*^^ be the set of multiaffine polynomials in 7. 

Theorem 2.1 (Part of Theorem 3.5 in PTj). Let T : ]R[x]^^'^ — )■ M[x] he a linear operator 
acting on the variables y:.. If the polynomial T{(x + yY"-^) G 6K[x,y], thenT maps ©]r[x]^^ 
into 6r[x]. 

Once the multivariate VT-Eulerian polynomials are shown to be stable, we can then reduce 
them to real stable univariate polynomials using the following operations. 

Lemma 2.2 (Part of Lemma 2.4 in [ST]). Given i,j G [n], the following operations preserve 
real stability of f E M[x].- 

(1) Differentiation: / i-^ |f-- 

(2) Diagonalization: / 1— )■ /|^.=^. ,. 

(3) Specialization: for a eM, f ^^ f\xi=a- 



Finally, there is an easy-to-clieck condition for real stability that we will use to show that 
certain polynomials are not real stable. 

Theorem 2.3 (Theorem 5.6 in [4\). Let f E M.[x]^^. Then f is real stable if and only if for 
all z,j e [n] and for all a G M", §-{a)£-{a) - ^(a)/(a) > 0. 

We note that most of these results have a complex counterpart, but for our purposes real 
stability suffices — all the polynomials we consider have positive integer coefficients. For this 
reason we will sometimes refer to real stable polynomials simply as stable polynomials. 

3. Stable ly-EuLERiAN Polynomials 

For clarity and completeness, we begin with a proof of the stability of the multivariate 
Eulerian polynomial of type A due to Branden. 

3.1. Eulerian Polynomials of Type A. Let An denote the Coxeter group of type A of 
rank n. We can regard A^ as Sym(?7, + 1), the group of all permutations on [n + 1] with 
generators S = {si, . . . , s„}, where Sj is the transposition (i, z + 1) for 1 < z < n. 

Proposition 3.1 (Proposition 1.5.3 in [2J). Given a = ai . . . cr„+i G An, 

^a(o-) = {si e S : ai > (Ti+i}, 

The following theorem is well-known. It was already known to Frobenius and has been 
reproven many times since. 

Theorem 3.2 (p. 829 of [12]). 

(1) An{x) = J2 ^"'"^'^^' 

cr&An 

has only real roots. 

Theorem 13.21 can be proved using RoUe's Theorem or by the method of interlacing roots, 
but the result becomes transparent by introducing additional variables and using properties 
of stability. 

Definition 3.1. Given a E An, define the type A descent top set to be 

T»T^(cr) = {max(o-j,o-i+i) : 1 <i <n,ai> o-j+i}, 

and similarly, let the type A ascent top set be 

AJAicr) = {max(cri,cri+i) ■.l<i<n,ai< ai+i}. 

For example, when a = 31452 G A4, I'T^(a) = {3,5} and ^T^(a) = {4,5}. Note the 
seemingly superfluous notation max(o"j, ctj+i) simply reduces to ctj and ctj+i in the case of 
type A descent top and ascent top sets, respectively. Its significance will become apparent 
when we introduce the type B descent top and ascent top sets. 

Theorem 3.3 (Branden [5]). 

(2) A„(x,y)=5^x^^-('^)y^^-W 

is stable. 



Proof. We proceed by induction. Note that Ao{xi,yi) = 1 is stable. By observing the effect 
of inserting n + 1 into a permutation a G An-i on the type A ascent top and descent top 
sets, we obtain the following recursion. For n > 0, we have 

(3) ^„(x, y) = {xn+i + ?/„+i)A„_i(x, y) + x„+i?/„+i9A„_i(x, y). 

We remind the reader here that d = XlILi {'^ ~^ W~) • ^^ ^^ easy to check using The- 
orem 12.11 that the linear operator T = (x„+i + Vn+i) + Xn+iVn+id is stability-preserving, 
because 

^1 1 A / 1 1 



r((x+u)W(y+v)[^ 



^n+lUn+l 



--E 



Vn+l -^n+l ■_-, \Xi ~r Ui 



(x+u)W(y+v)[ 



In J{_ when x,y,u,v£j{+ 



is in (5ik[x, y, u, v]. The result follows. D 

Specializing the yi variables to 1, it follows that 
Corollary 3.4. 

A„(x) = J2 ^^^^^"^ 

is stable. 

Diagonalizing x, we obtain 
Corollary 3.5. 

is stable. 

Since An{x) is univariate, this corollary is equivalent to the statement that An{x) has only 
real roots (Theorem 13. 2p . 

We refer to [15j for a proof of the stability of a (slightly different) multivariate refinement 
of the classical Eulerian polynomials. That refinement has close connections with the afiine 
Eulerian polynomial of type C, which we will address in more detail in Section 14. 2[ 

Next, we present our results. We start by defining the multivariate Eulerian polynomial 
of type B and proving that it is stable. 

3.2. Eulerian Polynomials of Type B. Let Bn denote the Coxeter group of type B of rank 
n. We regard -B„ as the group of all signed permutations on [±n] = {— n, . . . , — 1, 1, . . . , n} 
with generators S = {so,Si,..., s„_i}, where Sq is the transposition (—1,1) and Sj = {i,i + l) 
for 1 < i < n — 1. Given a = (cti, . . . , cr„,) G -B„, we let 

Nia) = \{ie[n]:ai<0}\ 

denote the number of negative entries in the signed permutation a. 

Type B descents have a simple combinatorial description that we will exploit. 
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Proposition 3.6 (Corollary 3.2 of [7], also Proposition 8.1.2 of [2]). Given a E Bn, 

l^sicr) = {si E S : ai > aj+i}, 

where ctq = 0. 

Analogously to type A, the type B Eulerian polynomials have only real roots. 
Theorem 3.7 (Brenti [7]). 
(4) Bnix) = J2 ^'^''^"^' 

has only real roots. 

In [7], Brenti introduced a "g-analog" of the univariate Eulerian polynomials and showed 
the following. 



TVH^IDsHI^ 



Theorem 3.8 (Corollary 3.7 of [7J). Forq>0, 

(5) Br,{x; q) = ^ q 

has only real roots. 

These Bn{x] q) polynomials specialize to Eulerian polynomials An-i{x) and Bn{x) — when 
g = and Q' = 1, respectively — so that Theorem 13.81 simultaneously generalizes Theorems 13.21 
and l3.7[ 

We will proceed in the same way that Theorem 13.31 extends Theorem 13. 2[ Recall that the 
stability of the multivariate refinement of the type A Eulerian polynomials in (|2]) came from 
the choice of the statistic. Choosing the larger index from each ascent and descent allowed for 
the simple stability-preserving recursion. We apply this idea to signed permutations in such 
a way that the definitions remain consistent with the definitions for ordinary permutations. 

Definition 3.2. Given a E Bn, define the type B descent top set to be 
DTij(a-) = {max(|cri|, \(Ti+i\) : < z < n - 1, a^ > ctj+i}. 

Analogously, we define the type B ascent top set to be 

AlBicr) = {max(|cri|, |cri+i|) : < i < n - 1, ai < ctj+i}. 
For example, when a = (3, 1, -4, -5, 2) E B^, T'TB(a) = {3, 4, 5} and ^TB(a) = {3, 5}. 

Theorem 3.9. For q> 0, 

(6) S„(x, y;q)=Yl g^('^)x'^^^('^)y^^«(-) 

is stable. 

Proof. As in the proof of Theorem 13. 3^ we proceed by induction. Bi{xi, yi]q) = qxi + yi is 
stable when q > 0, which settles the base case. By observing the effect on the ascent top and 
descent top sets of type B of inserting n + 1 or — (n + 1) into a signed permutation a E -B„, 
we obtain the following recursion. For ri > 0, we have 

(7) S„+i(x, y; q) = {qXn+i + i/„+i)S„(x, y; g) + (1 + g)a;„+i?/„+i9fi„(x, y; q). 
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To complete the proof, we note that for a fixed q > 0, the hnear operator acting on the 
right hand side, T = (gx„ + ?/„) + (1 + q)xnynd preserves stabihty by Theorem 12 ■![ since 



r((x+u)W(y+v)[^ 



^n+l?/n+l 



1 v^ / 1 + ^ 1 + g 



E 



(x+u)'»I(y+v)'" 



is in ©r[x, y, u, v] whenever g > 0. D 

Our theorem has some noteworthy consequences. The value of i?„(x, y;g) at g = — 1 is 
immediate from the recursion. 

Corollary 3.10. 

S„,(x,y;-l) = (y-x)W 

If we set g = 1, we obtain the analogue of Theorem 13.31 for type B. 



Corollary 3.11. 



is stable. 



i?„(x,y;l)=5^x^^-My^^-('^) 



(tG-B„ 



We would like to point out that when we plug in g = into ([H]), we get a homogenized 
polynomial that is not equal to the polynomial 74„_i(x, y) from ([2]), since their recursions 
differ. Rather, i?„(x, y;0) is the permanent of the following n x n matrix M = [rriij) 
considered in [6]. For i,j G [n], let rriij = Xi, when i < j and rriij = yj, otherwise. When we 
expand the permanent by the last column, we obtain the recurrence in ([7]) with q = (see 
Lemma 3.3 in [6] for a proof). 

Specializing the y variables in 5„(x, y; q) to 1, it follows that 
Corollary 3.12. For q > 0, 

5.(x;5)= j:g^(-)x^^-(-) 

is stable. 

Finally, observe that (the non- homogeneous) -Bn(x; q) does reduce to 74„_i(x) and -Bn(x) — 
the multivariate Eulerian polynomial of type A and type B — when q = and g = 1, re- 
spectively. Diagonalizing x in Bn{^]q) yields the polynomial Bn{x]q) defined in ([5]). We 
therefore recover Theorem 13.81 as a corollary. 

Corollary 3.13. For q > 0, Bn{x;q) is stable. 

3.3. Eulerian Polynomials for Colored Permutations. Theorem 13.91 can be extended 
in two directions simultaneously: from signed permutations to colored permutations, and 
from a single q variable to several. 

Let Tjr denote the cyclic group of order r with generator (. We will take ( to be an rth 
primitive root of unity. The wreath product G^ = ZrlAn-i is the semidirect product (Z.^)^"" xi 
An_i. Its elements can be thought of as a = (C^^^i, • • • , C'^"''n), where Cj G {0, 1, . . . , r — 1} 
and r G A„_i. G^ is sometimes called the generalized symmetric group. Its elements are 
also known as r-colored permutations, which reduce to signed permutations and ordinary 
permutations when r = 2 and r = 1, respectively. In other words, -B„ = G^ and A„_i = G^. 
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Definition 3.3. Given a = {C'^^Ti, . . . , C^"Tn) G G^, let J^I"(o") be the multiset in which each 
i G [n] appears Cj times. 

Note that for a G Bn, |^(cr)| = N{a), the number of negative entries in cr = (ai, . . . , an)- 
We adopt the following total order on the elements of (Z^ x [n])U {0} (see [IJ, for example): 

C~^n <---<Cn< C~\n-1) <■■■< C{n-1) <■■■ < C'H <---<Cl<0<l<---<n. 

Using this ordering, the definitions of descent, descent top set, and ascent top set all 
extend verbatim from Bn to G^. We shall therefore, by slightly abusing the notation, use 
the same symbols to denote them. For example, when a = (3, C^l, C^4, ("^5, (2) G G^, then 
DTB(a) = {3,5}, andyiTij(a) = {3,4,5}. 

Branden generalized Brenti's Bn{x;q) polynomial, defined in ([5]), to multiple q variables, 
and proved the following. 

Theorem 3.14 (Corollary 6.5 in [3]). Let q = (gi, ...,g„). If Qi > 0, for 1 <i <n, then 



)^|I'sWI_ 



(8) 5„(x;q)= 5^q^(- 

creBn 

has only simple real roots. 

Next, we extend this result simultaneously to G^^ and to multiple x variables. 

Theorem 3.15. If qi > 0, for all 1 < i < n, then the multivariate (\-Eulerian polynomial 
for the generalized symmetric group, G^^, defined as 

(9) G;(x,y;q)=5^q^(-)x^^-(-)y^--M 

is stable. 

Proof. G^{xi, yi] qi) = {qi + ■ ■ ■ + qi~^)xi + yi is clearly stable when gi > 0. The theorem 
follows immediately from the following recursion. For n > 1, 

Gli?^, y; q) = [((g„ + ■ ■ ■ + (i^~^)xn + 1/„) + (1 + ■ • ■ + g;-^)x„i/„a] G';_i(x, y; q). D 
As a consequence, we obtain a generalization of Corollary 13. 101 to G^^. 
Corollary 3.16. Let r > 2. For an rth root of unity, ( ^ 1, we have 

G;(x,y;£_^) = (y-x)W. 

n 

Letting r = 2 also generalizes Theorem 13.91 to multiple q variables. 

Corollary 3.17. If qi > 0, for all 1 < i < n, then i?„(x, y; q) is stable. 

Diagonalizing q gives us a result for G^ with a single q variable. 

Corollary 3.18. Ifq>0, then ^^(x, y;g) is stable. 
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4. Stable W^-Eulerian Polynomials 

Dilks, Petersen and Stembridge studied Eulerian-like polynomials associated to affine Weyl 
groups. They defined the so-called "affine" VT-Eulerian polynomials as the "affine descent" - 
generating polynomials over the corresponding finite Weyl group. In [9], they showed that 
the (univariate) W^-Eulerian polynomials have only real roots for types A and C, and also for 
the exceptional types. We strengthen these results for types A and C by giving multivariate 
stable refinements of these polynomials as well. 

4.1. AfRne Eulerian Polynomials of Type A. Let An denote the Coxeter group of type 
A of rank n. The affine descents of type A contain the (ordinary) descents of type A and 
an extra "affine" descent at if and only if, an+i > (Ji, where a = (a"i, . . . ,an+i) G An- 
Formally, 

T)A{a) = DA(a) U {sq : a„+i > ai}. 
See Section 5.1 in [9J for further details. 

The definitions of descent top and ascent top sets for type A can be extended in the 
obvious way. For a E An, 

V7A{(r) = DTA(a) U {a„,+i : a„+i > aj, 

A7A{cr) = A7A{(r) U {ai : a^+i < ai} 
and we obtain the following result. 

Theorem 4.1. _ _ 

I„(x,y)=5^x^^-(-)y^^-(-) 

is stable. 

Proof. This statement is immediate once we establish the following lemma. Stability follows, 
since A„(x, y) is stable and the operator on the right-hand side is clearly stability-preserving. 

D 

Lemma 4.2. For n > 0, we have 

y4„(x, y) = (n + l)a;„,+i|/„+iA„_i(x,y). 

Proof. Consider a permutation a G An-i. We will modify it to obtain a permutation in An. 
Append (n + 1) to the end of a and pick a cyclic rotation of the newly obtained permutation. 
The new permutation will have the same affine ascent top and affine descent top sets as the 
ascent top and descent top sets a had and in addition it will have {n + 1) both as an ascent 
top and as a descent top. To conclude the proof, note that there are exactly n + 1 cyclic 
rotations. This is essentially a refinement of the proof of Proposition 1.1 of [16]. D 

By diagonalizing x and specializing y to 1 , Lemma 14.21 reduces to an identity discovered 
by Fulman (Corollary 1 in [13]) and we recover the result of Dilks, Petersen, and Stembridge: 

Corollary 4.3 (see Section 4 of [H]). 



Mx) = Yl 



X 

has only real roots. 
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One can construct a recurrence from Lemma 14.21 for the afiine Eulerian polynomials of 
type A as well, but this recurrence will not preserve stability. 

4.2. AfRne Eulerian Polynomials of Type C. Let C„ denote the Coxeter group of type 
C of rank n. Affine descents of type C consist of the ordinary descent set of type C, which 
coincides with the descent set of type B (see Proposition 13.61 for type B descents) and an 
extra "afiine" descent at when cr„ > 0. Formally, 

Vcia) = Vcia) U {so : a„ > 0}. 

See Section 5.2 of [9] for further details. 

As in type A, the definition of the type C afiine ascent and descent top sets can be adapted 
from those of type C (equivalently, type B). For a & Cn, let 

V7cia) = VJBicr) U {a„. : a„ > 0}, 
AJcicr) = A7Bi(T) U {an : a„ < 0}. 

Theorem 4.4. _ _ 

C„(x,y)=5^x^^-My^^-(-) 

is stable. 

Proof. Ci{xi, yi) = 2xiyi and the following recurrence holds for n > 1: 

C„(x,y) = 2x„y„5C„_i(x,y). D 

We note that a very similar recurrence also appeared in [15| (without the factor of two and 
with a different initial value), in connection with stable Eulerian polynomials over Stirling 
permutations. 

There is also a direct connection between the polynomials C(x, y) and A(x, y). 

Proposition 4.5. 

Cn{Xi, . . . , Xn, yi, . ■ . , yn) = 2 XnynAn~l[XQ, . . . , X^-i, y^, . . . , ^n-l)- 

Proof. Follows by a similar argument as Lemma 14. 2[ D 

Once again, this gives a refinement of the univariate identity by Fulman [Tl] . We mention 
that multivariate refinements (different from the above) and respective refinements for the 
identities for A(x, y) and C(x, y) have appeared in 0. 

5. Towards Stable Refinements of D{x), B{x), D{x) 

5.1. Eulerian Polynomials of Type D. Let Dn denote the Coxeter group of type D of 
rank n. Recall that Dn can be thought of as the (order 2) subgroup of Bn consisting of 
all the signed permutations with an even number of negative entries. Specifically, Dn has 
generators S = {sq, Si, • • • , s^}, where Sq = (—2, 1)(— 1, 2) and Sj = (z,i + 1) for 1 < i < n. 

Proposition 5.1 (Proposition 8.22 of [2J). Given a E Dn, 

1)d{o-) = {si e S : a^ > ctj+i}, 



, def 

where ctq = — (T2. 
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Based on the real-rootedness results for types A, B, and the exceptional types (and some 
computer evidence) Brenti conjectured the following. 

Conjecture 5.2 (Conjecture 5.1 in [7j). The type D Eulerian polynomial has only real roots. 

This was — and still is — the remaining unproven part of Conjecture 11.11 Unfortunately, 
extending the type B ascent top and descent top definitions in the naive way does not work. 
If we let 



(10) 



'D7]j{(t) = {max(|crj|, |o"i+i|) :0<z<?t. — 1,(Tj> ctj+i}, and 
A7d{o-) = {max(|a-i|, |cr,+i|) : < i < n - 1, ai < ctj+i}, 



then the type D Eulerian polynomial is not multiafiine (for example, the monomial corre- 
sponding to cr = 123 is ylvz)- Furthermore, it fails to be stable. For n = 3, we have 

Dl (x, y) = xlx2. + 2x2X3y2 + x^yl + xly^ + Ax2X2.y2. + 4x3^/3 + 2x22/22/3 + 4x3^/22/3 + vlvz + ^x^vl- 

When ^2 = 2/3 = X3 = 2 + i, and X2 = (-1 + 2i){2i + y/3) ^ -5.7321 + 1.4641^ (which are all 
in the upper half plane), D3(x, y) = 0. (We use the D*(x, y) notation for the multivariate 
Eulerian polynomial of type D to avoid confusion with a different multivariate generalization 
which will appear later on.) 

Branden used the refinement given in ([8]) to prove an intriguing result about enumerating 
type B descents over type D permutations. 



Theorem 5.3 (Corollary 6.10 in [3]). 



has only real roots. 






X 



25bWI 



It would be interesting to find a multivariate analog of this result. Again, the straight- 
forward application of our method with the descent top set definition in flTU]) results in a 
non-stable polynomial already for n = 3. 

Table [T] gives the type B and type D descents for permutations in D3. As demonstrated 
above, our choice oi X = x\ao\ and Y = 2/|cro| does not result in stable polynomials, though it 
is possible that some other choice will. 



a 


x'^'^s 


^Wn 


a 


^Ws 


^Wn 


a 


^m. 


^m. 


a 


^■i^'J, 


^VJn 


123 


Yy2y3 


^2/22/3 


123 


Xx2y3 


-'^ 2/22/3 


123 


Xy2X3 


Yy2X3 


123 


Yy2X3 


Xy2X3 


132 


yysxs 


Yy3X3 


132 


Xxsys 


Xxsys 


132 


Xysxs 


Yy3X3 


132 


YxsVs 


Xx3y3 


213 


Yx2y3 


Yx2y3 


213 


^2/22/3 


^2/22/3 


213 


X2/2X3 


Xy2X3 


213 


YX2X3 


YX2X3 


231 


Yysxs 


Yysxs 


231 


Xx^ys 


Xxsys 


231 


Xy3X3 


Yy3X3 


231 


Yx3y3 


Xx3y3 


312 


Yx3y2 


Yx3y2 


312 


^2/32/2 


Xy-iy2 


312 


Xy3X2 


Yy3X2 


312 


YX3X2 


Xx3a;2 


321 


YX3X2 


YX3X2 


321 


^2/32/2 


^^2/32/2 


321 


Xy3X2 


X2/3X2 


321 


Yx3y2 


Yx3y2 



Table 1. Type B and type D descents over the permutations in D^ 
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Another result related to resolving Conjecture 15.21 came from [8], in which Chow found 
a recurrence for Dn{x). Sadly, the resulting formula is much more complicated than its 
counterparts in types A and B, and does not seem amenable to a multivariate generalization. 
We reproduce the recurrence here, fixing a typo from the original paper (the term in the box 
was mistakenly written with a minus sign). 

Theorem 5.4 (cf. Theorem 5.3 in [8J). 

Dn+2ix) = (n(l + 5x)+4s)D„+i(x) 
+4x(l -x)D'^^-^{x) 

+ ((1 - xf - n{l + 3xf - 4:n{n - 1)2;(1 + 2x))Dn{x) 
-{Anx{l - x)(l + 3x) + 4x(l - xf)D'^{x) 
-4x\l-xyD';{x) 



+ {2n{n - 1)2;(3 + 2x + 3x^) +4n(n - l)(n - 2)x^(l + x) )Dn-i{x) 

+ {2nx{l - xf{3 + x) + 8n{n - 1)2;2(1 - x){l + x))D'^_^{x) 
+4^2(1 - x)2(l + x)D'^_^{x). 

Encouraged by the simplicity of our methods for types A and B, we propose a new line of 
attack. Stembridge showed that the Eulerian polynomials of types A, B, and D are related 
via the following identity. 

Theorem 5.5 (Lemma 9.1 of [19j). For n > 2, 

(11) Dn{x) = Bn{x) - n2"-ixA„__2(a;). 

By replacing the univariate polynomials Bn{x) and An{x) in ( ITT1) by their stable multi- 
variate generalizations given in (jj]) and ([6]), respectively, we obtain the following multivariate 
refinement of Dn{x): 

(12) ^„(x,y) = S„,(x,y;l) -n2"-ix„y„A„_2(x,y). 
For example, when n = 2, we obtain 

D2{xi,X2,yi,y2) = B2{xi,X2,yi,y2;l) - 4x22/2 ^o(2;i,?/i) = {xi + 2/i)(a;2 + 2/2), 
and when n = 3, the polynomial is 

1)3 (x, y) = 2:1X2X3 + X2X32/1 + XiX3y2 + X3yiy2 + X1X22/3 + a^22/i2/3 + a:i2/22/3 + 2/i2/22/3 
+4(x2X32/2 + xixsys + X22/22/3 + xsyiys). 

These polynomials fail to be stable, even for n = 3. This follows from Theorem 12.31 since 
specializing the y variables to 1 gives 

D3(x) = 1 + Xi + X1X2 + X1X2X3 + 5(X2 + X3 + X1X3 + X2X3) 

and 

^^3(x) ai^3(x) a^i^3(x) ^ , . .„ 

which fails to be nonnegative for X2 > 0. It is possible that a different refinement would 
result in a family of stable polynomials, from which real-rootedness would follow. 
It does seem to be the case, however, that the coefficients remain non-negative. 
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Conjecture 5.6. D.„(x, y) G N[x, y] for all n>2. 

We have verified this conjecture by computer for n < 11. If true, then (fT2|) suggests a 
refinement of the descent statistic of type D which in turn could lead to a new recursion for 
the type D Eulerian polynomials. 

5.2. AfRne Eulerian Polynomials of Types B and D. Dilks, Petersen, and Stembridge 
noted that the only missing cases from Conjecture ll.2l are types B and D. They considered 
multivariate refinements of these polynomials, but they indexed the variables by descents 
and not by descent tops. It was noted in [15] that polynomials indexed by the descents fail 
to be stable for type A. 

Two new identities relating ordinary and affine Eulerian polynomials also appeared in [9]. 

Proposition 5.7 (Proposition 6.1 in [9j). For n > 2, 

2Cn{x) = Bn{x) + 2ntCn-i{x). 

This identity could be used to obtain a multivariate refinement for the unsettled type B 
case in a similar way as suggested for type D above. Let 

(13) S„(x, y) = 2C'„(x, y) - 2nx„2/„5„_i(x, y; 1) . 

In contrast to the type D case, these polynomials turn out to be stable for n < 4 (we have 
verified this by computer calculations). Therefore, we suggest to investigate whether this is 
a stable multivariate of the afiine Eulerian polynomial of type B. 

Conjecture 5.8. 5„(x,y) G ©K[x,y], for n > 2. 

Also, this multivariate refinement seems to be monomial positive. 

Conjecture 5.9. Bn{:s., y) G N[x, y], for n > 2. 

Finally, for sake of completeness, we give another identity by Dilks, Petersen and Stem- 
bridge which relates the following polynomials. 

Proposition 5.10 (Proposition 6.2 in [9J). For n > 3, 

(14) Bn{x) = Dn{x) + 2ntD„_i(z). 

This identity might be helpful for finding a multivariate refinement of the affine Eulerian 
polynomial of type D. 

6. Conclusion and further remarks on the statistics A7, DT 

In this paper, we have extended the stability results for the multivariate type A Eulerian 
polynomial to Eulerian polynomials (and affine Eulerian polynomials) of some other Coxeter 
groups and the generalized permutation group Gn,r- A cruical step in our proofs was to find 
a suitable generalization of the descent top and ascent top statistics to these groups. 

The descent statistic {i : a^ > o"i+i} for permutations has a long history, going back to 
the works of Carlitz and Riordan. Recall that for a permutation of length n there are 2"~^ 
possible descent (and ascent) sets, as each position i E [n — 1] can either be a descent or an 
ascent. In particular, for any subset S C [n — 1] there is a permutation whose descent set is 
S. Finding the number of permutations with a given descent (and ascent) set is a classical 
problem — one uses a standard inclusion-exclusion argument (see Example 2.2.4 in |17j). 
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We now find the possible descent and ascent top sets (of type A), wliicli corresponds to 
tlie possible monomials in A„(x, y). (From here on, I)'J and A7 refers to the descent top 
and ascent top sets of type A.) 

Proposition 6.1. Gtven {1)7, AT) C [n] x [n] define ? = T)7nA7 andV = [n] \(DTU>AT). 
The pair (2)T, A7) is the descent top and ascent top set of some permutation of [n] if and 
only if, for all i = I, . . . ,n, \[i]r\V\> \ [i] n 7\ . 

Proof. Clear from the recurrence given in ([3]). D 

In fact, the {'D7,A7) pair of statistics turns out to be equivalent to a well-known statistic 
studied by Frangon and Viennot [llj, called "type" — the quadruple of statistics consisting 
of the set of peaks, valleys, double descents, and double ascents. In particular, the set 
of peaks {ctj : crj_i < Uj > ctj+i} is CP, defined above, and similarly, the set of valleys 
{(Tj : CTj-i > (Jj < CTi+i} is V. For the set of double descents 2)D := {cr,j : ctj-i > cr,j > ctj+i} 
we have that 2)2) = 'D7\A7 (and similarly for the double ascents we have DA = A7\'D7). 

In light of the above, the following theorem should not be a surprise. 

Theorem 6.2. The number of distinct monomials in the expansion o/v4„_i(x, y) is counted 
by Cn, the nth Catalan number. 

Proof. We give a bijection from (2)T, A7) G [n] x [ri\ to 2-colored Motzkin paths of length 
n — 1. A 2-colored Motzkin path of length n — 1 is a lattice path of N^ running from (0, 0) to 
{n — 1, 0) that never goes below the x-axis and whose allowed steps are the NE step (1, 1), 
the SE step (1, —1) and the colored E step (1, 0) which is labeled by one of two colors. For 
i = 2, . . . ,n, let the (j — l)st step of the Motzkin path be SE, NE, E, or E_ depending 
on which set (IP, V, 2)2), or DA, respectively) j belongs to. The 2-colored Motzkin paths of 
length n — 1 are known to be counted by the Catalan numbers C„ (see, for example, [18j). D 

Exploring this connection, we can also give a formula for the coefficient on a particular 
monomial x^'^y^'^ — that is, the number of permutations with a given descent and ascent top 
set. Following Viennot, we define a valuation on these 2-colored Motzkin paths. NE and 
SE steps starting at height A; > are given a weight a^ = c^ = fc + 1. We give E steps a 
weight of 6fc = k + 1 and E_ steps a weight of 6^ = k + 1, which is equivalent to forgetting that 
we had two different types of E steps and weighting an (uncolored) east step by bk = 2k + 2. 
Viennot's theory of Laguerre histories [20], which gives a combinatorial interpretation for 
the moments of the Laguerre polynomials, allows us to then recover the permutations with 
given descent and ascent top sets. In particular, the number of permutations associated to 
a path is simply the product of the weights of the steps in the path. 

Similar reasoning works for type B as well. 

Theorem 6.3. The number of distinct monomials in the expansion of Bn{^,y; 1) is counted 
by Cn, the nth Catalan number. 

Likewise, the coefficents on a particular monomial in i?„(x,y; 1) are given by Viennot's 
theory by using the weights a^ = 2{k + 1), 6^ = 4A; + 2, and Ck = 2k. It is a simple extension 
to show that a monomial x^^y^^ in -B„(x, y; q) has a g-coefficient given by weighting NE 
steps by a^ = (A; + 1)(1 + q), E steps by 6a: = 1 + k{l + q), E_ steps by 6^ = g + k{l + q), and 
SE steps by c^ = k{\ + q). It is also easy to extend these results to G'„,r(x, y; q). 
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Appendix A. List of stable multivariate Eulerian polynomials 

^o(x,y) = l 

Ai(x,y) =X2+y2 

A2(x, y) = X2X3 + X3t/2 + a^2l/3 + 2X3?/3 + |/2?/3 

Si(x,y) =xi +1/1 

-B2(x, y) = xi^s + X2yi + xii/2 + 4a;2?/2 + l/i?/2 

-Bi(x,y;g) = ?a;i + ?/i 

S2(x, y; g) = g^XiX2 + qx2yi + qxiy2 + (1 + qfx2y2 + l/i?/2 

Ai(x,y) =2x22/2 

^2(x, y) = 2x2X-iy-i + 2x^y2y-i 

A3(X, y) = 2X2X3X4?/4 + 2X3X4l/2?/4 + 2X2X42/31/4 + 4X3X42/3?/4 + 2X42/2?/32/4 

Ci(x,y) =2xi2/i 

C2(X, y) = 4xiX2l/2 + 4X2?/l|/2 

C'3(x, y) = 8xiX2X3?/3 + 8x2X3^12/3 + 8x1X32/21/3 + 16x2X32/22/3 + 8x32/12/22/3 
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